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The  modelling  effort  involved  development  of  user-oriented  software  that  is 
compatible  with  USATHAMA  computers  for  prediction  of  heatup  rates  and  agent 
migration  and  degradation  as  a  function  of  process  variables  with  a  selection 
of  boundary  conditions  to  represent  different  structural  conditions.  The 
concurrent  experimental  work  involved  (a)  development  of  a  method  for 
laboratory  simulation  of  HOP  and  monitoring  desorption  products  using  a 
time-of-f light  mass  spectrometer  with  a  molecular  beam  inlet  (TOFMS);  (b) 
spiking  and  analysis  of  concrete  with  HD,  GB,  and  VX  by  gas  chromatography;  and 
(c)  use  of  the  TOFMS  to  monitor  thermal  desorption  of  agent  spiked  concrete. 

The  model  developed  on  this  program  provides  a  rapid  and  convenient  means  for 
examining  the  thermal  consequences  of  variations  in  HDP  parameters. 

Analyses  of  concrete  core  samples  obtained  from  Newport  Army  Ammunition  Plant 
(NAAP)  showed  no  evidence  for  residual  VX.  Spiking  experiments  showed  that  all 
three  agents  were  degraded  by  contact  with  concrete,  but  at  different  rates  and 
with  different  consequences.  Mustard,  HD,  is  degraded  slowly  by  both 
hydrolysis  to  form  oxathiane  and  by  rearrangement  to  form  dithiane.  The  latter 
process  is  accelerated  thermally.  GB  reacts  rapidly  with  the  basic  components 
of  concrete  by  replacement  of  the  fluorine  group  with  a  hydroxyl  group.  VX  is 
degraded  to  a  complex  spectrum  of  products.  This  process  is  initially  slow, 
but  appears  to  be  autocatalytic,  i.e.  at  least  some  of  the  products  accelerate 
the  degradation  process.  In  all  three  cases,  many  of  the  degradation  products 
are  relatively  volatile  and  are  readily  desorbed  from  the  concrete  matrix. 

The  current  work  has  verified  some  of  the  degradation  chemistry  proposed  by 
other  studies  of  agents,  and  suggests  that  the  degree  of  contamination  of 
concrete  structures  may  depend  not  only  on  the  extent  of  the  original 
contamination  but  also  on  the  history  of  aging  of  the  structure.  In  cases  of 
relatively  casual  original  contamination,  the  long  term  effects  of  the  chemical 
reactions  with  the  concrete  may  have  essentially  decontaminated  the  structure. 
This  was  apparently  the  case  with  the  core  samples  obtained  from  NAAP,  where  no 
residual  VX  was  detected.  On  the  other  hand,  experience  with  severely 
contaminated  HD  sites  has  shown  that  appreciable  residual  HD  can  remain  trapped 
with  the  structure  after  many  years.  The  slowness  of  the  HD  degradation 
coupled  with  possible  saturation  of  the  reaction  capacity  of  the  concrete  could 
account  for  this  observation.  It  is  therefore  recommended  that  consideration 
of  application  of  HDP  to  structural  decontamination  be  preceded  by  a  formalized 
core  sampling  and  analysis  program  to  aid  in  customizing  the  application  of  the 
process  to  indi'^dual  cases.  It  is  also  recommended  that  the  model  be  expanded 
to  include  cher  :al  factors  developed  on  this  program  to  further  facilitate 
customization  of  the  process. 
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1.0  Introduction 


As  part  of  their  responsibilities  in  DOD  real  property  disposal,  the 
U.S.  Army  Toxic  and  Hazardous  Materials  Agency  (USATHAMA)  must  develop  and 

demonstrate  technology  to  identify,  contain,  and  eliminate  toxic  and  hazardous 
materials  and  related  contamination.  In  support  of  this  mission,  Battel le's 
Columbus  Division  (BCD)  has  recently  completed  two  projects  that  identified 
promising  decontamination  systems  on  a  bench  scale  in  the  laboratory.  The 

objective  of  the  current  project  is  to  further  define  the  fate  of  chemical 

agents  (HD,  VX,  and  GB)  adsorbed  on  structural  components  such  as  concrete  and 
paint  films.  The  project  goals  are  to  elucidate  the  chemistry  and  kinetics  of 
such  interactions  through  (a)  development  of  a  kinetic  model  for  desorption  and 
transport  of  agents  from  porous  substrates  and  paints;  (b)  conduct  of 

experiments  with  agent  simulants  to  verify  the  kinetic  model;  and  (c)  conduct 
of  experiments  with  agents  to  demonstrate  the  utility  of  the  model  and  to 
identify  specific  chemical  species  desorbed  and  residual  products  of  the 
process. 

The  overall  approach  included: 

1.  Development  of  a  user-oriented  transient  heat  and  mass  flow 
model  for  use  with  80286/80287  based  computers  with  provision 
for  selection  between  multiple  boundary  conditions  choices  and 
multiple  output  modes. 

2.  Use  of  simulants  (half-mustard  and  di isopropylmethylphosphonate) 
doped  on  concrete  to  develop  methodology  for  study  of  thermal 
desorption  of  agents  from  concrete  using  a  time-of-f light  mass 
spectrometer  with  a  molecular  beam  sampling  inlet  (TOFMS) . 

3.  Analysis  of  core  samples  obtained  from  NAAP. 

4.  Spiking  and  gas  chromatographic  (GC)  analyses  of  concrete 
samples  doped  with  HD,  GB,  and  VX. 
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2.0  PROCEDURES 


2.1  MODELING  OF  HDP 


The  hot  gas  desorption  process  involves  exposing  one  or  both  sides  of 
concrete  wall  to  a  flowing  stream  of  gas  at  elevated  temperatures  in  order  to 
heat  the  wall  to  temperatures  sufficient  to  cause  decomposition  and/or 
desorption  of  adsorbed  agent.  As  depicted  in  Figure  1,  practical  application 
of  the  process  can  involve  a  variety  of  heat  transfer  processes,  including 
transport  to  the  surface  of  the  waP  (or  to  both  sides  of  the  wall  in  some 
instances);  unsteady  conduction  within  the  wall;  and  heat  losses  through  the 
wall  to  either  the  outside  air  or  to  the  soil.  This  can  further  be  complicated 
by  the  presence  of  an  insulating  barrier  between  the  wall  and  the  air.  Mass 
flow  within  the  system  is  similarly  complicated  by  the  fact  that  degradation, 
evaporation,  and  condensation  can  all  occur  within  the  heated  wall. 

The  purpose  of  this  portion  of  the  overall  task  was  to  develop  a  user 
oriented  computer  code  that  is  sufficiently  flexible  to  deal  with  the  multiple 
variations  of  the  process  and  which  will  allow  ready  coupling  to  prediction  of 
residual  agent  concentrations.  The  following  discussion  outlines  the  general 
approach  taken  to  the  solution  of  the  heat  and  mass  flow  problem,  and  discusses 
use  of  the  code.  A  detailed  description  of  the  mathematics  of  the  problem  is 
given  in  Appendix  A. 

2.1.1  Computational  Approach 

The  most  convenient  approach  to  maintain  flexibility  in  choice  of 
process  variables  and  boundary  conditions  is  through  the  use  of  a  finite 
element  description  of  the  process.  Several  different  procedures  are  available 
for  this  approach,  as  described  by  Lienhardt  (1981),  Crank  (1956),  Carnahan 
(1969),  and  many  other  standard  texts  on  heat  and  mass  transfer.  These 
procedures  all  involve  solution  of  a  series  of  simultaneous  equations  that  are 
generated  by  expressing  the  fundamental  differential  equations  in  terms  of 
Taylor  series  expansions  of  the  differentials.  This  results  in  solutions  at  a 
finite  intervals  (nodes)  within  a  temperature/depth  or  concentration/depth 
array.  It  is  generally  agreed  that  the  Crank-Nicolson  method  (for  detailed 
discussion  see  Crank,  op.cit.)  yields  the  best  accuracy  per  step  size  used  in 
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the  calculation.  For  example,  the  accuracy  of  the  Crank-Nicolson  method  is 
proportional  to  the  square  of  the  number  of  spatial  nodes  used,  whereas  in  the 
implicit  method  used  in  the  Battel! e  Phase  II  report  (1985),  the  accuracy  is 
directly  proportional  to  the  number  of  spatial  m  '<  Furthermore,  fundamental 
mathematical  considerations  of  the  Crank-Nicolson  method  indicate  that  it  is 
appropriate  to  use  time  steps  that  are  directly  proportional  to  the  square  of 
the  spatial  step,  a  fact  that  greatly  simplifies  the  equations  used. 

2.1.2  Code  Description 

The  code  for  the  HDP  model  was  written  in  C  (Microsoft  v5.0)  for  use 
with  80286/80287  processors.  It  is  menu-driven  with  user  selected  (1)  gas 
properties,  (2)  structural  materials  and  dimensions,  (  boundary  conditions, 
(4)  contaminant  type,  (5)  contaminant  distribution  anc  vel ,  d  (6)  output 
modes.  As  examples,  the  user  can  select  any  of  four  bom  ary  iditions  (both 
sides  of  wall  heated,  one  side  air  cooled,  one  side  insu  or  one  side  in 
contact  with  moist  soil);  and  output  can  be  tabular  or  <e  form  of  either 
time  or  depth  profiles  of  temperature,  adsorbed  phas  ncentration,  or  gas 
phase  concentration.  The  code  supports  hard  copy  outpu.  co  plotters  ing  HPGL 
or  bit-image  graphics  to  either  IBM  or  Epson  compatible  printers.  ile  a 
default  set  of  recommended  materials  properties  is  included  in  the  code, 
virtually  every  parameter  can  be  changed  to  simulate  special  cases. 

In  its  current  stage  of  development,  only  the  heat  transfer  portion 
of  the  HDP  model  is  functional.  Previsions  have  been  made  for  the  mass 
transfer  modules  and  peripheral  I/O  segments  related  to  the  mass  transfer  have 
been  incorporated,  but  time  and  funding  limitations  related  to  the  need  for 
incorporation  of  key  data  from  the  experimental  portion  of  the  program  have 
dictated  postponement  of  completion  and  testing  of  the  full  program.  The 
accuracy  of  the  heat  transfer  module  (temporarily  de-coupled  from  mass 
transfer)  has  been  checked  extensively  with  each  boundary  condition  and  a 
variety  of  changes  in  input  parameters.  In  the  case  of  both  sides  of  the  wall 
being  heated,  the  model  results  were  compared  with  the  explicit  solution  that 
is  available.  In  the  cases  of  the  other  boundary  conditions,  calculations  were 
compared  with  other  machine  solutions  that  were  available  for  specific  sets  of 
conditions.  Extensive  evaluations  also  were  made  by  both  Battelle  and  USATHAMA 
personnel  for  ease  of  use  and  general  utility  of  the  code,  and  recommendations 
for  improvement  of  useability  have  been  incorporated  into  the  overall  code 
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design.  The  code  also  has  been  optimized  for  speed  with  the  most  complex  case 
(determined  by  the  64k  per  data  segment  limitation  of  the  80286  processor) 
being  completed  in  about  40  sec  (Compaq  286). 

Text  introducing  the  use  of  the  coo^  is  included  on  the  diskette 
containing  the  executable  image  file  of  the  code. 

2.2  CORE  ANALYSIS  AND  SPIKIN'!  EXPERIMENTS 


2.2.1  Newport  Core  Analyses 

A  set  of  concrete  samples  was  received  from  the  NAAP  for  analyses  of 
VX  and  VX  degradation  products.  The  core  removed  from  the  trench  area  was 
suspected  of  having  the  greatest  degree  of  contamination.  The  other  cores  were 
also  taken  from  areas  where  potential  exposure  could  have  occurred  except  for 
the  background  core  removed  from  the  loading  area  of  the  facility. 

Each  of  the  cores  was  divided  into  approximately  1  inch  segments,  and 
the  first  (top  surface)  and  third  segments  were  pulverized  using  a  diamond  mill 
to  a  finely  divided  powder.  A  small  sample  of  the  powder  from  each  segment 
(10  +  0.1  g)  was  weighed  into  a  4-dram  vial  tor  extraction.  The  samples  were 
extracted  using  chloroform  (10  +  0.1  ml),  and  vortex  mixing  was  utilized  for 
30+5  seconds  to  facilitate  the  contact  between  the  concrete  particles  and  the 
solvent.  The  solution  was  allowed  to  settle,  and  the  liquid  layer  was 
transferred  through  a  0.45-/«n  syringe  filter  to  a  gas  chromatograph  (GC)  vial. 
The  samples  were  analyzed  using  capillary  GC  with  a  flame  photometric  detector 
(GC/FPD).  The  analog  signal  from  the  detector  was  integrated  using  Beckman 
Peak  Pro  Chromatography  software. 

2.2.2  Agent  Spiking  Experiments 

The  agent  spiking  experiments  were  conducted  to  determine  if  the 
agents  of  interest,  sarin  (GB) ,  mustard  (HD),  and  VX,  were  degraded  with  time 
when  exposed  to  concrete.  These  experiments  were  conducted  by  adding  agent  to 
1-g  quantities  of  pulverized  concrete  from  the  NAAP  core  samples.  The  concrete 
cores  were  pulverized  using  a  diamond  mill  which  produces  a  finely  divided 
powder.  Two  types  of  experiments  were  conducted:  dilute  agent  experiments  and 
concentrated  agent  experiments. 
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In  the  dilute  agent  experiments,  initially  dilute  agent  solutions  of 
agent  in  hexane  were  prepared.  The  concentration  of  these  solutions  was  2  mg 
of  agent/ml  of  solution  for  each  agent.  One  gram  c  the  pulverized  concrete 
core  samples  was  added  to  5  ml  of  each  dilute  agen.  solution.  The  solutions 
were  then  mixed  by  swirling  and  allowed  to  stand  until  the  concrete  particles 
had  settled  to  the  bottom  of  the  solution.  At  this  time,  an  aliquot  of  the 
liquid  layer  was  removed  from  each  solution  and  analyzed  by  6C/FPD  to  confirm 
the  agent  concentration.  Subsequent  aliquots  were  removed  at  multiple  time 
intervals  to  determine  the  status  of  the  agent  concentration.  If  the  agent 
concentration  was  observed  below  the  measurable  concentration  (0.02  percent), 
the  solution  was  fortified  with  10  mg  of  neat  agent  and  rean  lyzed. 

In  the  concentrated  agent  experiments,  1  ml  of  neat  chemical  surety 
material  (CSM)  was  added  to  1  gram  of  pulverized  concrete  core  sample.  Again, 
aliquots  of  the  liquid  were  removed  at  multiple  time  intervals  and  analyzed  by 
GC/FPD  to  determine  the  agent  concentration.  A  gas  chromatograph-mass 
spectrometry  (GC-MS)  analysis  also  was  conducted  on  concentrated  samples  after 
the  samples  had  aged  in  an  attempt  to  identify  the  degradation  products. 


2.3  THERMAL  DESORPTION 


2.3.1  Apparatus 

A  schematic  representation  of  the  TOFMS  used  for  the  thermal 
desorption  studies  is  shown  in  Figure  2.  This  unit  consists  of  a  CVC  Model  MA- 
3A  time-of-f light  mass  spectrometer  that  is  equipped  with  a  molecular  beam 
sampling  inlet  and  a  sample  chamber.  The  sample  chamber  is  attached  to  the 
molecular  beam  sampler  via  an  0-ring  flange.  Within  the  sample  chamber,  a 
heater  and  sample  cup  assembly  is  mounted  on  a  centrally  located  moveable 
pedestal.  A  thermocouple  us^d  for  measuring  sample  temperature  feeds  through 
the  pedestal  to  the  base  of  the  cup.  In  use,  the  height  of  the  sample  cup  is 
adjusted  to  1-2  mm  below  the  sampling  orifice.  Air  drawn  through  one  of  the 
ports  on  the  sample  chamber  at  approximately  30  cc/min  sweeps  over  the  sample 
and  into  the  orifice.  The  first  stage  of  the  sampler  is  maintained  at  a 
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pressure  of  approximately  0.1  torr.  Sample  gas  entering  this  region  undergoes 
a  rapid  expansion,  with  only  those  molecules  traveling  along  a  straight  line 
between  the  sampling  orifice  and  the  skimmer  orifice  entering  the  mass 
spectrometer.  The  ion  source  of  the  mass  spectrometer  also  lies  along  this 
straight  path  so  that  the  sampling  process  is  essentially  collisionless  and 
preserves  the  integrity  of  the  species  being  sampled.  Calibration  of  the 
sensitivity  of  the  apparatus  is  performe  I  by  placing  samples  of  neat  liquids 
having  known  vapor  pressures  in  the  sample  cup.  The  principal  components  of 
the  sample  air  (nitrogen,  oxygen,  argon,  etc.)  serve  as  mass  markers  for 
calibration  of  the  mass  range.  Because  of  the  dependence  of  the  sampling  rate 
through  the  orifice  on  temperature,  the  sensitivity  decreases  slightly  with 
temperature,  but  a  typical  detection  limit  for  DIMP  is  0.042  /* g/cc. 

2.3.2  Method  Development 

Prior  to  performing  the  agent  desorption  experiments,  preliminary 
experiments  were  conducted  with  the  simulants  half-mustard,  HM,  and 
di isopropylmethylphosphonate,  DIMP.  These  preliminary  experiments  involved 
spiking  1  cm  X  1  cm  X  0.6  cm  high  porosity  (17  voi  K)  and  low  porosity 
(12  vol  %)  concrete  coupons  with  the  simulants.  The  spiking  solutions  used 
consisted  of  HM  and  DIMP  mixed  with  hexane  at  a  concentration  of  60  mg  of 
simula.it/ml  of  solution. 

Each  spiked  sample  was  loaded  onto  the  sample  cup  in  the  sample 
chamber.  The  sample  chamber  was  then  attached  to  the  molecular-beam  inlet  of 
the  TOF.  A  background  spectrum  was  obtained  just  before  mounting  the  sample 
chamber.  The  temperature  of  the  sample  was  then  programmed  from  ambient  to 
200-250°C  at  rates  of  2-10  degrees/min. ,  and  the  mass  range  40  to  200  was 
repeatedly  scanned.  In  analyzing  the  mass  scans,  typically  the  two  masses  that 
were  used  as  calibration  masses  were  40  (Argon)  and  one  of  the  major  masses  of 
the  simulant. 

2.3.3  Agent  Spiked  Samples 

All  of  the  agent  experiments  were  conducted  using  1  cm  X  1  cm  X  0.6 
cm  high  porosity  (17  vol  %)  concrete  coupons.  These  coupons  came  from  the  same 
high  porosity  concrete  batch  that  was  used  in  the  simulant  experiments.  The 
concrete  samples  were  spiked  with  one  of  three  agents:  GB,  HD,  or  VX.  The 
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purity  of  each  agent  was  checked  before  spiking.  The  purities  of  the  GB,  HO, 
and  VX  were  found  to  be  96  percent,  97  percent,  and  84  percent  respectively, 
which  are  within  normal  limits.  Both  fresh  and  aged  samples  were  prepared  for 
each  agent.  The  fresh  samples  were  prepared  by  applying  25  /il  of  the  agent  to 
top  surface  of  a  concrete  coupon  immediately  before  it  was  analyzed  with  the 
mass  spectrometer. 

Prior  to  preparing  the  aged  samples,  a  concrete  coupon  was  soaked  in 
water  for  24  hours  to  determine  the  approximate  amount  of  water  a  coupon  would 
absorb.  By  weighing  the  coupon  before  and  after  soaking,  it  was  determined 
that  the  coupon  absorbed  7  percent  of  its  weight  in  water.  It  was  decided  that 
50  percent  of  the  absorbed  water  weight  would  be  added  to  the  aged  coupons; 
therefore  water  equivalent  to  3.5  percent  of  the  sample  weight  was  added  to 
each  aged  coupon.  After  the  addition  of  water,  100  /tl  of  agent  was  applied  to 
each  coupon.  Coupons  spiked  with  the  same  agent  were  then  placed  in  a  closed 
container  and  allowed  to  age  prior  to  the  sample  being  analyzed  with  the  mass 
spectrometer. 

Before  analyzing  an  agent  spiked  sample,  a  calibration  run  was 
performed  with  the  agent  to  determine  the  sensitivity  of  the  mass  spectrometer 
for  the  agent  and  to  identify  a  typical  mass  scan  for  the  agent.  Here  a  small 
stainless  steel  plate  was  placed  in  the  sample  cup.  Next  a  small  amount  of  the 
agent  was  placed  on  the  plate  and  the  sample  chamber  was  mounted  on  the  mass 
spectrometer.  Mass  scans  were  taken  at  room  temperature  and  at  elevated 
temperatures.  From  the  mass  scan  for  the  agent  itself,  the  following 
sensitivity  factor  was  calculated: 

S  =  (AMT  Agent) *(IAr)  (1) 

(AMT  Ar)*(Iagent) 

where, 

S  *  Sensitivity  factor 

AMT  Agent  *  partial  pressure  of  agent  at  room  temperature 
AMT  Ar  =  partial  pressure  of  argon  at  room  temperature 
IAr  »  Intensity  of  most  prominent  Ar  peak 
Iagent  =  Intensity  of  most  prominent  agent  peak 

Following  the  calibration  run,  samples  spiked  with  the  calibration 
agent  were  analyzed.  Before  and  after  the  desorption  run  on  a  spiked  sample,  a 
background  scan  was  obtained  with  the  sample  chamber  removed.  In  a  typical 
desorption  run,  the  sample  was  loaded  onto  the  sample  cup,  and  the  sample 

10 


chamber  was  mounted  on  the  molecular-beam  inlet  of  the  TOF.  The  temperature  of 
the  sample  was  then  increased  at  2-10  degrees/min.  from  near  ambient  to 
200-277°C  or  until  no  further  desorption  occurred.  On  some  occasions  during 
heating,  samples  were  held  at  one  temperature  for  several  minutes  to  study 
certain  desorption  products  before  heating  was  begun  again.  While  the  GB  and 
HD  spiked  samples  were  heated,  the  mass  range  of  40  to  200  was  repeatedly 
scanned.  For  the  VX  spiked  samples,  the  mass  range  of  40  to  280  was  repeatedly 
scanned.  In  analyzing  the  mass  scans  of  the  agent  spiked  samples,  the  masses 
used  for  calibration  masses  were  40  (Argon)  and  149  (a  moderate  intensity 
background  peak) . 

After  performing  the  desorption  runs  an  attempt  was  made  to  estimate 
the  total  amount  of  agent  recovered  from  the  sample.  This  estimate  was  done 
using  the  equation  below: 

N  =  (AMT  Agent)* (At)  =  S*(Iagent*At)*(AMT  Ar)*M 

(IAr)  (2) 

N  =  Amount  of  agent  recovered,  grams 

(Iagent*At)  =  Area  under  the  peak  intensity  vs  time  plot  for  agent 


M 


=  Molecular  weight  of  agent 


3.0  RESULTS 


3.1  HOP  MODEL 


Figure  3  shows  the  main  menu  of  the  HOP  code;  several  sub-menus  are 
used  depending  on  selections  made  from  the  main  menu.  Figure  4  shows  two 
examples  of  the  hard-copy  output  capabilities  of  the  HDP  code. 

3.2  CORE  ANALYSIS  AND  SPIKING  EXPERIMENTS 


3.2.1  Newport  Core  Analyses 

Essentially,  all  core  sample  analyses  were  identical  with  no  evidence 
of  VX  or  its  degradation  products  being  observed  above  the  1  /*g/g  detection 
limit  of  the  GC/FPD  detector.  The  continuation  of  these  analyses  was  stopped 
since  the  most  likely  candidates  for  contamination  were  found  free  of  any 
measurable  concentrations  of  VX  or  chromatographable  degradation  products. 

3.2.2  Agent  Spiking  Experiments 

A  summary  of  the  GC/FPD  analyses  for  the  dilute  agent  experiments  is 
presented  in  Table  1.  Figure  5  shows  a  graph  of  the  results  from  the  dilute 
experiment  conducted  on  HD.  The  half-life  of  the  HD  has  been  indicated  on 
Figure  5.  As  can  be  seen  on  Table  1,  HD  was  found  to  be  much  more  stable  than 
VX  or  GB  in  the  dilute  agent  experiments.  Both  VX  and  GB  showed  evidence  of 
rapid  decomposition  when  mixed  with  the  pulverized  concrete. 

Because  the  HD  appeared  to  be  relatively  stable  in  the  dilute 
experiments,  only  VX  and  GB  were  subjected  to  the  concentrated  agent 
experiments  in  which  1  ml  of  neat  agent  was  added  to  1  gram  of  pulverized 
concrete.  The  results  of  the  GC/FPD  analyses  from  the  concentrated  experiments 
are  shown  in  Table  2.  These  results  show  that  even  undiluted  VX  and  GB  neat 
agent  decomposed  rapidly  when  applied  to  the  pulverized  concrete.  In 
comparison,  the  undiluted  GB  sample  decomposed  more  rapidly  than  the  undiluted 
VX  sample  because  virtually  no  GB  was  remaining  after  168  hours  but  all  of  the 
VX  was  still  remaining  after  264  hours.  Figure  6  shows  a  graph  of  the  results 
from  the  concentrated  VX  experiment.  Along  the  100  percent  plateau  shown  in 
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Figure  3.  Main  Menu  of  HDP  Code 


Hot  Gas  Properties 

Gas  (F)low  Rate  =  5000  cfm  (G)as  Temperature  =  1000.0  F 

Wall  Properties 

(C) onductivity  =  0.70  Btu/ft-hr-F  (S)pecific  Heat  =  0.21  Btu/lb-F 

(D) ensity  *  144  1 b/f t3  (T)hickness  =  8.0  in 

(H) eight,  Width  *  12.0  x  12.0  ft  (I)nitial  Temperature  *  70.0  F 

Boundary  Conditions 

(I) .  Both  sides  heated  (2).  Air  cooled  side 

(3).  Moist  soil  side  (4).  Insulated  side 

Output  Choices 

(5).  IBM  Printer  (6).  Plotter  (7).  Disk 

Chemical  Agent 

(8).  HD  (9).  GB  (0).  VX 


Set  up  problem  by  pressing  letter  or  number  key  for  each  variable  you 
wish  to  change  or  set. 
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Depth,  i  riches 

4A.  Hard-Copy  Output  of  Temperature  as  a  Function 
Depth  at  Several  Different  Times  in  a  Concrete 


Figure  4B.  Hard-Copy  Output  of  Temperature  as  a  Function  of 

Time  at  Several  Different  Depths  in  a  Concrete  Wall 


TABLE  2.  SUMMARY  OF  ANALYSIS  RESULTS  FOR  THE 
CONCENTRATED  NEAT  AGENT  CONCRETE  STUDIES 


(1  ml  CSM  on  1  gram  Concrete) 


Date  of  Sampling  Agent 

Hours  After 
Initial  Dose 

Percent 

Remaining 

09/01/88 

VX 

0 

100 

09/14/88 

VX 

264 

100 

10/07/88 

VX 

816 

0.3 

11/02/88 

VX 

1368 

0.2 

11/08/88 

VX 

1536 

0.3 

11/16/88 

GB 

0 

100 

11/23/88 

GB 

168 

0.04 

12/19/88 

G8 

792 

0.0* 

*Be1ow  detectable  limit  (0.02  percent). 
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Figure  6,  the  normally  clear  VX  liquid  changed  to  a  yellow  colored  liquid. 
Along  the  0  percent  plateau  on  Figure  6,  the  liquid  remaining  was  found  to  be 
viscous.  A  GC-MS  analyses  was  conducted  on  the  concentrated  VX  sample  after 
the  sample  had  aged  1824  hours  in  an  attempt  to  identify  any  degradation 
products  that  were  present.  The  results  of  the  analysis  are  shown  in  Appendix 
B. 


3.3  THERMAL  DESORPTION  EXPERIMENTS 


3.3.1  Simulant  Desorption  Runs 

Table  3  contains  a  summary  of  the  simulant  desorption  runs.  As  noted 
in  Table  3,  the  HM  samples  that  were  treated  with  loadings  of  600  ug/cm2  HM  (in 
hexane)  and  allowed  to  dry  yielded  no  detectable  desorption  products.  The  only 
HM  spiked  concrete  samples  that  yielded  desorption  products  were  those  that 
were  either  heavily  spiked  and  were  analyzed  immediately  without  allowing 
hexane  to  evaporate  completely,  or  those  that  were  actually  soaked  in  the  HM 
solution  overnight.  Even  in  the  last  sample  in  Table  3,  where  the  concrete  was 
exposed  to  HM  vapor  in  a  closed  jar,  there  was  no  evidence  for  retention  of 
either  HM  or  its  reaction  products.  Rather,  the  residual  HM  was  found  only  in 
the  air  contained  in  the  jar. 

With  DIMP,  the  re  ults  were  similar  to  the  HM  results  in  that  the 
DIMP  was  observed  to  desorb  only  with  freshly  prepared  samples.  No  desorption 
products  besides  the  DIMP  itself  were  identified  on  the  DIMP  spiked  samples. 
The  freshly  prepared  sample  yielded  a  desorption  curve  as  a  function  of 
temperature  as  shown  in  Figure  7.  The  slope  of  this  curve  represents  Ep/R,  the 
apparent  energy  of  the  process  divided  by  the  universal  gas  constant.  By 
multiplying  the  slope  by  R  (.001987  kcal/mole-K)  it  was  found  that  the  apparent 
energy  of  the  process  Ep  was  16.5  +  1.1  kcal/mole.  Measurements  of  the  vapor 
pressure  of  DIMP  yields  the  following  Antoine  equation  for  the  vapor  pressure 
of  DIMP  (courtesy  of  E.  C.  Penski,  CRDEC) : 

(43.5-102  C)  Log io  P  -  5.7202  - 878 • 624  — 

(110.578  +  t,  C)  (3) 

Several  vapor  pressures  for  DIMP  were  calculated  using  Equation  3. 
These  pressures  were  then  plotted  with  respect  to  temperature  in  Figure  8.  The 
slope  of  the  curve  in  Figure  8  represents  AHV/R  (AHV  =  heat  of  vaporization  of 
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TABLE  3.  SUMMARY  OF  SIMULANT  DESORPTION  RUNS 
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HM/Hexane  Vapor  None  detected 

Vapor  plus 
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DIMP) .  By  multiplying  the  slope  by  R,  a  AHy  of  17.8  kcal/mole  was  found  for 
DIMP.  The  similarity  of  the  values  for  Ep  and  AHV  indicate  that  the  desorption 
of  DIMP  from  the  concrete  basically  involved  the  evaporation  of  the  DIMP. 

With  the  use  of  Equation  2,  an  estimate  was  made  on  the  amount  of 
DIMP  that  was  recovered  from  the  freshly  prepared  sample.  In  particular,  it 
was  found  that  only  approximately  6  percent  of  the  DIMP  was  recovered.  This 
low  recovery  amount  indicates  that  the  DIMP  was  evaporating  very  quickly  and 
that  most  of  the  DIMP  had  evaporated  during  the  short  time  span  between  when 
the  sample  was  spiked  and  when  the  sample  chamber  was  loaded  onto  the 
molecular-beam  sampling  head. 

Overall,  the  results  indicated  that  the  samples  spiked  at  the  surface 
with  HM  and  DIMP  did  not  retain  the  simulants  for  any  appreciable  period  of 
time,  and  that  the  loss  of  the  spike  can  occur  by  reaction  and  or  evaporation 
at  room  temperature.  These  results  indicated  that  excessive  loading  or 
exposures  over  extended  periods  of  time  would  be  required  in  the  agent 
desorption  experiments,  especially  for  HD,  in  order  to  achieve  significant 
penetration  into  the  concrete. 

3.3.2  Results  with  Agents 

Table  4  shows  a  summary  of  the  agent  desorption  runs.  Figures  9-15 
show  peak  intensity  and  temperature  plotted  against  time  for  all  of  the 
desorption  runs. 

3.3.2. 1  GB  Runs 

For  GB,  the  only  desorption  product  found  on  the  sample  freshly 
spiked  with  25  ul  of  agent  was  GB  itself.  For  the  aged  samples  in  which  water 
had  been  added,  the  major  desorption  product  was  found  to  be  isopropyl  methyl 
phosphonate  (IMP)  with  a  small  amount  GB  still  bring  present.  This  result 
indicates  that  both  time  and  moisture  are  required  ior  the  GB  to  react  with  the 
concrete.  In  the  two  aged  samples,  the  IMP  was  first  detected  in  the 
temperature  range  26-45°C  which  is  within  the  range  of  ambient  temperatures  to 
which  a  concrete  structure  could  be  exposed.  For  the  GB  sample  aged  for  120 
hours,  GB  was  first  detected  at  98°C,  while  for  the  GB  sample  aged  for  336 
hours,  GB  was  first  detected  at  181°C.  The  increased  temperature  required  to 
detect  GB  in  the  older  sample  suggests  that  the  depth  of  the  reaction  layer  has 
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TABLE  4.  SUMMARY  OF  AGENT  DESORPTION  RUNS 
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Time  Cmin.) 

Figure  9.  Plot  of  Peak  Intensities  ( m/e  =  99) 
for  Fresh  GB  Sample. 


GB  on  Concrete  (336  hrs) 
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Time  (min) 

Figure  11.  Plot  of  Peak  Intensities  (m/e  =  97  and  99) 
for  336  Hour  GB  Sample. 


HD  on  Concrete  (freshly  treated) 
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Figure  12.  Plot  of  Peak  Intensities  (m/e  =  109)  for 
Fresh  HD  Sample. 
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Time  (min) 

Figure  13A.  Plot  of  Peak  Intensities  (m/e  =  104,  109  and  120) 
for  312  Hour  HD  Sample. 
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Figure  14A.  Plot  of  Peak  Intensities  (m/e  =  104,  109  and  120) 
for  360  Hour  HD  Sample. 
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increased  in  the  older  sample  so  that  the  6B  has  penetrated  deeper  into  the 
concrete  matrix. 

The  formation  of  the  IMP  probably  results  from  the  GB  reacting  with 
hydroxide  (OH-1)  ions  present  in  the  moist  concrete.  This  reaction  is  shown 
be  I  w  in  Equation  4.  The  hydroxide  ion  (OH-1)  appears  to  simply  substitute  for 
the  fluoride  ion  (F-l)  in  the  GB  structure.  The  mass  spectra  for  GB  and  IMP 
were  found  to  be  identical  with  the  exception  that  the  IMP  spectrum  was  shifted 
down  two  mass  units.  This  shift  can  be  accounted  for  by  the  substitution  of  an 
OH-1  for  the  F-l  in  GB.  Any  major  restructuring  of  the  molecule  would  be 
evidenced  by  more  significant  differences  in  the  two  mass  spectra.  Figure  16 
shows  a  comparison  of  the  mass  spectra  for  the  GB  and  IMP,  while  Figure  17 
shows  a  typical  mass  spectrum  from  the  desorption  of  the  aged  samples  in  which 
both  GB  and  IMP  are  present. 

o  a 

II  8 

H  C  —  P  —  OCHCCH^D  +  OH  >  H  C  —  P  —  OCHC CH  D  +  F 

3  |  3  2  ^  3  |  3  2 

F  OH 

(4) 

GB  IMP 

Using  equation  2  an  estimate  was  made  of  the  amount  of  GB  recovered 
from  the  three  GB  spiked  samples.  Approximately  85  percent  of  the  25  ul  of  GB 
applied  to  the  freshly  spiked  sample  was  recovered.  This  estimate  seems 
reasonable  because  some  GB  was  lost  between  the  time  the  sample  was  spiked  ana 
the  sample  chamber  was  loaded  onto  the  molecular-beam  inlet.  For  the  two  aged 
samples,  recovery  estimates  were  made  by  considering  the  sum  of  the  GB  and  IMP 
peaks.  Fifty-seven  percent  of  the  GB  appeared  to  be  recovered  from  the  120 
hour  GB  sample,  and  30  percent  of  the  GB  appeared  to  be  recovered  from  the  336 
hour  GB  sample. 

3. 3. 2. 2  HD  Runs 

As  was  found  with  GB,  the  only  desorption  product  identified  for  the 
freshly  spiked  HD  sample  was  the  agent  itself.  For  the  two  aged  HD  samples, 
the  desorption  products  found  were  HD,  1,4-dithiane,  and  1,4-oxathiane.  Of 
these  three  desorption  products,  HD  was  by  far  the  most  abundant  product  as  can 
be  seen  on  Figures  13  and  14.  Between  the  1,4-dithiane  and  the  1,4-oxathiane, 
the  1,4-oxathiane  appeared  to  be  more  abundant.  Again  the  presence  of 
desorption  products  besides  that  of  the  agent  itself  on  the  aged  samples 
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Relative  Intensity  Relative  1,. tensity 


Figure  16,  Comparison  of  Mass  Spectra  for  GB  and  IMP 


Mass  to  Charge  Ratio 


indicates  the  importance  of  time  and  moisture  in  order  for  the  agent  to  react 
with  the  concrete.  Figure  18  shows  a  mass  spectrum  for  an  aged  HD  sample  in 
which  HD,  1,4-oxathiane,  and  1,4-dithiane  are  present. 

As  is  shown  on  Table  4,  HD  was  first  detected  in  the  temperature 
range  of  29-47°C  which  is  within  the  range  of  ambient  temperatures  to  which  the 
concrete  could  be  exposed.  On  the  aged  samples,  the  1,4-oxathiane  was  first 
detected  at  around  170°C  in  both  cases.  The  1,4-dithiane  was  first  detected  at 
the  same  temperature  as  1,4  oxathiane  in  the  312  hour  sample  and  at  a 
temperature  41°C  higher  in  the  360  hour  sample.  The  greater  temperature 
required  to  detect  the  1,4-dithiane  in  the  360  hour  sample  sugges's  that  the 
1,4-dithiane  may  have  been  located  deeper  in  the  concrete  matrix  on  the  360 
hour  sample  than  on  the  312  hour  sample. 

The  1,4-oxathiane  and  1,4-dithiane  are  cyclic  compounds  that  appear 
to  be  the  reaction  products  of  the  de-chlorination  reactions  which  occur  in  the 
concrete.  Shown  below  in  Figure  19  is  a  comparison  of  the  structure  of  HD  with 
the  ring  structures  of  1,4-oxathiane  and  1,4-dithiane. 


Figure  19.  Structural  Comparison  of  HD  and  its  Desorption  Products. 

The  formation  of  1,4-oxathiane  requires  hydrolysis  of  the  HD 
molecule.  Here  the  HD  reacts  with  hydroxide  (OH-1)  ions  and  eventually 
liberates  HC1  to  form  the  1,4-oxathiane.  Equation  5  shows  the  hydrolysis 
reaction. 
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The  formation  of  the  1,4-dithiane  results  from  the  thermal 
degradation  of  HD.  This  process  involves  the  association  of  2  HD  molecules 
with  the  elimination  of  dichloroethane  to  form  the  thermal  degradation  product 

1.4- dithiane.  Thermal  degradation  of  HD  to  form  1,4-dithiane  occurs  at  room 
temperature  but  is  accelerated  by  heating. 

From  equation  2,  it  was  estimated  that  85  percent  of  the  25  /*1  of  HD 
applied  to  the  freshly  spiked  sample  was  recovered.  As  with  the  freshly  spiked 
6B  sample,  this  estimate  seems  reasonable  because  some  HD  was  lost  between  the 
time  the  sample  was  spiked  and  the  sample  chamber  was  loaded  onto  the 
molecular-beam  inlet.  For  the  two  aged  samples,  recovery  estimates  were 
considered  using  only  the  HD  peaks.  No  attempt  was  made  to  use  the 

1.4- oxathiane  or  1,4-dithiane  peaks.  On  the  two  aged  samples,  39  percent  HD 
recovery  was  found  on  the  312  hour  sample,  and  50  percent  HD  recovery  was  found 
on  the  360  hour  sample. 

3. 3. 2. 3  VX  Runs 

The  analysis  of  desorption  products  from  VX  spiked  samples  was  found 
to  be  difficult.  In  particular,  once  the  VX  was  heated  above  140°C,  a  liquid 
condensed  on  the  skimmer  orifice  and  plugged  it  so  that  no  further  mass  scans 
could  be  taken.  An  analysis  of  this  liquid  with  the  mass  spectrometer 
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indicated  the  liquid  was  VX.  As  a  result  of  the  VX  condensation  problem,  only 
one  run  was  conducted  on  a  VX  spiked  sample.  Modifications  such  as  mounting  a 
heater  inside  the  molecular-beam  sampling  head  would  have  to  be  made  to  the 
present  apparatus  before  the  desorption  products  of  VX  could  be  analyzed  at 
temperatures  above  140°C. 

The  one  VW  sample  analyzed  was  a  sample  that  had  aged  48  hours.  On 
this  sample,  the  only  desorption  product  found  before  the  skimmer  orifice 
plugged  at  142°C  was  VX  itself.  As  can  be  seen  on  Figure  15,  no  significant  VX 
was  detected  until  the  sample  temperature  reached  87°C.  The  VX  peak  intensity 
continued  to  increase  at  temperatures  above  87°C  until  the  orifice  became 
plugged  at  142°C.  Figure  20  shows  a  typical  mass  spectrum  from  the  aged  VX 
sample.  No  estimates  were  made  of  the  amount  of  VW  recovered  because  of  the 
limited  number  of  points  that  could  be  measured. 


4.0  CONCLUSIONS AND  RECOMMENDATIONS 


The  HDP  model,  in  its  current  form,  provides  a  convenient  and  rapid 
means  for  examining  the  thermal  consequences  of  variation  in  HDP  parameters. 
However,  although  provisions  were  made  for  incorporation  of  chemical  parameter, 
the  model  ■%  currently  incomplete  with  respect  to  description  of  chemical 
release  rates.  It  is  recommended  that  consideration  be  given  to  completion  of 
the  HDP  model. 

Based  on  the  current  spiking  experiments,  it  is  concluded  that  all 
three  agents  (HD,  GB,  and  VX)  are  degraded  in  the  presence  of  concrete.  The 
degradation  of  HD  is  relatively  slow,  with  apparent  half-life  of  about  1747 
hours.  Indeed,  this  process  may  be  slowed  further  by  consumption  of  the  basic 
components  of  the  concrete.  The  principal  product  of  reaction  of  HD  with 
concrete  is  oxathiane.  Dithiane  also  was  observed  in  the  thermal  desorption 
experiments,  but  this  substance  is  believed  to  be  formed  by  a  normally  slow 
rearrangement  process  of  the  mustard  which  is  accelerated  by  the  heating 
process.  That  is,  the  dithiane  does  not  result  from  reaction  with  the 
concrete.  The  reaction  of  GB  with  concrete  is  rapid  and  involves  exchange  of 
the  flourine  group  of  the  GB  with  hydroxyl  groups  derived  from  the  basic 
components  of  the  concrete.  The  reaction  of  VX  with  concrete  is  slow 
initially,  but  the  process  of  degradation  of  VX  appears  to  be  autocatalytic.  A 
complex  mixture  of  products  is  formed  by  this  process.  The  current  work  has 
resulted  in  chromatographic  isolation  of  43  such  products.  GC/MS 
identification  of  these  products  is  currently  tentative  at  best  because  of  the 
paucity  of  agent  related  compounds  in  the  current  NBS  data  base.  It  is 
recommended  that  a  more  extensive  examination  of  the  samples  and  spectra  from 
the  degraded  VX  samples  be  performed  using  reference  spectra  available  from 
other  agencies  concerned  with  agent  analysis  and  using  standard  chemicals  for 
generation  of  reference  spectra  where  such  information  is  not  available. 
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APPENDIX  A 

rHE  MATHEMATICS  OF  HDP 


THE  MATHEMATICS  F  HDP 


The  HDP  involves  a  series  of  heat  and  mass  flow  elements: 

1.  heat  flow  to  and  mass  flow  from  a  hot  surface 

2.  heat  flow  coupled  with  vaporization/condensation  and 
pyrolysis  within  the  structure 

3.  heat  and  mass  flow  from  the  cool  side  of  the  structure 

The  core  of  the  problem  in  forming  a  mathematical  description  of  the 
process  lies  in  item  2,  whereas  the  boundary  conditions  are  derived  from 
consideration  of  items  1  and  3.  In  the  following  discussions,  abbreviated 
terms  are  explained  within  the  text  but  all  fundamental  terms  are  defined  in 
the  glossary  of  symbols. 


CONSERVATION  EQUATIONS 

At  any  point  within  the  wall,  conservation  of  heat  requires  that 


ks  ^  ^  -  qe  (  n>e  -  ii»a  )  -  qrkr  (Ca  +  Cg  ) 


(1) 


(total  heat)  =  (sensible  heat)  -  (heat  of  evap/cond)  -  (heat  of  pyrolysis) 


For  the  adsorbed  phase,  conservation  of  mass  requires 

(?) -I  (*-)(*-*■)  - 

(surface  concentration)  -  (evaporation/condensation)  -  (pyrolysis) 
For  the  vapor  phase,  conservation  of  mass  requires 


(s“)  '  D  (0s ) +  )  ("e'Aa )  "krC9  (3) 

(concentration)  *  f  diffusive  \  +  f  evaporation/  \  -  (pyrolysis) 

\  transport  /  \  condensation  / 

In  these  equations,  the  rates  of  evaporation,  me,  and  condensation,  ma,  are 
given  by 

ii»a  31  where  7  2  e-ASy/3R  (4) 


' T*TTT" "!■'!  1  .V ’tH' !'-") s*/v ;q •'  " 1,1  ^t1?^ 


and 


NsMv 


e-AHv/RT 


4RT 

The  rate  constant  for  pyrolysis,  kr,  is 
kr  -  Ae’E*/RT 

The  solution  to  the  problem  is  gained  by  simultaneous  solution  of 
equations  1-3  along  with  the  appropriate  initial  and  boundary  conditions. 


(5) 

(6) 


Boundary  Conditions 

To  provide  flexible  use  of  the  HOP  model,  4  different  boundary 
conditions  are  considered. 

1.  Both  sides  of  wall  heated. 

2.  One  side  of  wall  cooled  by  ambient  air. 

3.  One  side  in  contact  with  moist  soil. 

4.  One  side  covered  with  insulation. 

The  general  approach  used  requires  determining  "he  effect  of  each  of 
these  boundary  conditions  on  the  temperature  gradient  at  the  surface  of  the 
wall.  In  general,  we  seek  to  find  an  expression  of  the  form 

(£-).■ 

for  each  of  the  boundary  conditions. 

For  a  heated  boundary,  there  is  no  heat  loss  at  the  boundary,  so, 


ks  (£-).-  h(TR-Ts) 

where  h  is  the  overall  heat  transfer  coefficient  which  .  'omposed  of 
contributions  from  radiation,  natural  convection,  and  forctu  convection,  ie. 


h  s  hp  +  (h^nc  +  h2fc) 


1/2 


Similarly,  the  boundary  condition  for  an  air  cooled  side  is 


(£-).  =  h»(Ws) 


(9) 

(10) 
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When  one  side  is  insulated, 


*>(£-) 

V  dx  /  . 

where  k'  and  L'  are  the  conductivity  and  thickness  of  the  insulation 
respectively.  In  the  case  of  one  side  in  contact  with  moist  soil,  the  soil  is 
a  semi-infinite  medium,  ie.  the  depth  of  soil  affected  by  heat  transfer  from 
the  wall  increases  with  time.  Lienhardt  (ref)  suggests  representing  the  depth 
of  thermal  penetration  by 

Ld  -  (at) 1/2  (12) 

where  a  is  the  thermal  diffusivity  of  he  soil.  The  boundary  equation  then 
becomes 


k| 

L' 


(  T*  -  Ts  ) 


(11) 


kS 


(  Td  ’  TS  ) 


(13) 


Heat  and  Mass  Transfer  Parameters 

Forced  Convection: 

kaNu 
d 


hfc 


where 


ka  a  conductivity  of  hot  gas  *  1.845  x  10-4  (  ~  )  (  — -  ) 

v  pr  '  \  non  / 

cal,  cm3,  K“1 

Nu  =  Nusselt  number  »  0.023  Re®*^  Pr*^  (dimensionless) 


0.75 


298 


Re  3  Reynolds  number  =  =  --—ft— \  (dimensionless) 

d/*  dTg  \  Ig  / 

Q  3  volumetric  flow  rate,  cm3  sec-1 

d  =  hydraulic  diameter  of  heated  space  (taken  normal  to  gas  flow 
direction),  cm 

H  3  gas  viscosity,  g  cm-1  sec-1 
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1  ^jR  '  t'tt  ^  r  . ;  »* 


r,r^ '  ■* 


Pr  »  Prandl  number  a  - £fi -  (dimensionless) 

0.0855  +  Cp 

Cp  -  heat  capacity  of  gas  »  (6.449  =  1.4113  x  10-3Tg  -  8.07  x  10-8T2g)/29, 
cal  g-lK-1 

Tg  «*  gas  temperature,  K 

Natural  Convection: 

hmc  ■  0.525  ^  ^  (PrGr)1/4 

where  Gr  *  Grashof  number  *  Tg  -  Ts 

ft 

p  s  -  —  f  ^  ^  =  coefficient  of  volume  expansion  * 

p  \  dT  / 

7.51  x  10“3  -  l.i  x  10-5  Tg  +  6.63  x  10-9  j2g  .  1.4  x  10-12i3g,  k-1 
g  =  acceleration  of  gravity  *  980.7,  cm  sec2 

Radiation: 

hr  =  <re*  (l4g  -  t4s  )/  (  Tg  -  Ts  ) 

a  *  Stefan-Boltzmann  constant  *  1.355  x  10-12,  cal  sec-1  cm-2  k-4 
c*  =  effective  emissivity  =  0.37  (dimensionless) 


)  (dimensionless) 


Agent  and  Degradation  Product  Parameters 


The  following  table  summarizes  transport  and  degradation  parameters 
of  agents  and  known  degradation  products. 


SUMMARY  OF  CHEMICAL  PARAMETERS  RELEVANT  TO  HDP 


Chemical 

D@298k 

cm2/sec 

Vapor  Pressure 

■ 

Pyrolysis 

ASy 

eu 

AHy 

kcal/mole 

log  A 
sec~i 

E*  1 

kcal/mole  m 

HD 

29.4 

13.4 

9.25 

25  1 

GB 

.066 

27 

10.9 

8.18 

23.2  M 

VX 

.065 

29.3 

16.3 

8.98 

27.8  * 

DEMP 

.066 

28.8 

11.5 

1 

I 

B 

I 

1 

I 

I 

I 

B 


GLOSSARY  OF  SYMBOLS 


Symbol 

Meaning 

A 

Arrhenius  pre-exponential  factor 

a 

thermal  diffusivity  of  soil 

J> 

coefficient  of  expansion 

Ca 

concentration  of  adsorbed  phase 

Cg 

concentration  of  vapor  phase 

Cp 

heat  capacity  at  constant  pressure 

D 

effective  diffusivity  of  vapor 

d 

hydraulic  diameter 

E* 

activation  energy 

e 

porosity  of  concrete 

e* 

emissivity 

g 

acceleration  of  gravity 

7 

collisional  reaction  efficiency 

Gr 

Grashof  number 

h 

overall  heat  transfer  coefficient 

hfc 

forced  convection  heat  transfer  coefficient 

hnc 

natural  convection  heat  transfer  coefficient 

hr 

radiation  heat  transfer  coefficient 

AHV 

enthalpy  of  vaporization 

k' 

thermal  conductivity  of  insulation 

kr 

rate  constant  for  pyrolysis 

ks 

thermal  conductivity  of  concrete 

L 

thickness  of  concrete 

L' 

thickness  of  insulation 
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GLOSSARY  OF  SYMBOLS  (Continued) 


Ld 

ma 

me 

P 

Ns 

Nu 

Pr 

Q 

Qe 

qr 

r 

R 

Re 

P 

a 

ASV 

T 

t 

Ta 

Td 

t9 

Ts 


thickness  of  thermal  penetration  in  soil 
mass  rate  of  adsorption 
mass  rate  of  vaporization 
viscosity 

mole  fraction  of  sorbate  on  surface 

Nusselt  number 

Prandl  number 

volumetric  flow  rate 

heat  of  vaporization  per  unit  mass 

heat  of  pyrolysis 

average  pore  size  in  concrete 

gas  constant 

Reynolds  number 

density 

Boltzmann  constant 
entropy  of  vaporization 
absolute  temperature 
time 

air  temperature 
soil  temperature 
hot  gas  temperature 
surface  temperature 


v  mean  molecular  speed 

x  distance 
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APPENDIX  B 

GC-MS  ANALYSES  OF  CONCENTRATED  VX  ON  CONCRETE  SAMPLE 
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GC-MS  ANALYSES  OF  CONCENTRATED  VX  ON  CONCRETE  SAMPLE 


GC-MS  analyses  were  conducted  on  the  liquid  taken  from  the  sample  in 
which  1  ml  of  VX  was  applied  to  1  gram  of  pulverized  concrete.  These  analyses 
were  conducted  after  the  sample  had  aged  1824  hours  in  an  attempt  to  identify 
the  degradation  products  of  VX. 

Figure  1  shows  the  GC  scan  that  was  obtained  from  the  liquid.  The 
mass  spectra  that  were  identified  for  each  GC  peak  were  subjected  to  a  Hewlett- 
Packard  probability  matching  (PBM)  search  routine  in  an  effort  to  identify  the 
possible  degradation  products.  This  routine  searches  the  NBS  library  and  uses 
probability  matching  to  identify  the  compounds  whose  mass  spectra  best  match 
the  mass  spectra  found  for  each  GC  peak.  Table  1  shows  a  list  of  the  possible 
product  compounds  that  were  identified  for  each  of  the  19  major  GC  peaks.  It 
should  be  noted  that  the  compounds  shown  in  Table  1  only  provide  the  best  mass 
spectral  match.  However,  additional  testing  would  be  required  to  determine  if 
the  peak  retention  times  of  the  identified  compounds  match  the  peak  retention 
times  found  for  the  sample  liquid.  Only  by  matching  both  retention  time  and 
mass  spectra  can  a  positive  identification  of  the  degradation  products  be  made. 
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IC  of  DRTR: VXDEG2ND. 


Concrete  Sample  Aged  1824  Hours 


TABLE  B-l .  POSSIBLE  DEGRADATION  PRODUCTS  IDENTIFIED  IN 
PBM  MASS  SPECTRA  SEARCH 
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APPENDIX  C 

Source  Code 
for 

Hotgas  Desorption  Model 
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SOU,  F  CODE 


#include  <math.h> 
linclude  <float.h> 
linclude  <stdio.h> 

#include  <conio.h> 

linclude  <graph.h> 

linclude  <ctype.h> 

linclude  <malloc.h> 

linclude  <process.h> 

char  far  *buffer;char  buf 1 [10] ; 

int  BFL[10],i,jj,j,k,NP;float  pi [15] ; 

float  q®5000.0,*float  ks*. 0029; float  cp».21;float  tg®810.77; float  l*?r\32; 

float  ti*294.111;float  kp». 000207; int  ST=12;float  STZ; 

float  ts=285. 78; float  NX.ta.re.pr.nu.hfc.ka.gra.ha.kd.ad; 

float  he*12;float  wi*12;float  d;int  ypix=350; 

float  rho*2.3; float  rhod,cpg,cpd;  float  wa=25.0; 

float  far  th [43] [355]; 

float  yp,delx,delt,hy,hx,hr,gr,hc,lp,n,x,tt,dth,ht;int  nz; 
float  far  cs [43] [355] ; float  far  eg [43] [355]; 
char  resp[2];FILE  *stream;char  presp[2J;char  mode; 
char  fname[12];  char  res;char  buf[90];mt  far  tm[640];int  m[12]; 

Idefine  LPT1  0; 

Idefine  C0M1  0; 
struct  rccoord  rcoord; 
char  cres; float  coav; 
double  a,b,c; 


main() 


start:  for  (i=0;i<=9;i *•+)  BFL[i]=3; 

menu  (); 

matrix:  printf  ("%s" , "Number  of  hours  to  simulate  "); 
scant  &n)  * 

ST  »pow(350*rho*cp*pow(l  ,2)/(n*ks*3600) ,  .5)-2;if(ST>=42){ST=42;} 
printf  ("%s  %d", "Maximum  number  of  spatial  nodes  *  ",ST); 
a*modf ( ST*2 . 54/ 1 , &b ) ; ST* ( i nt ) ( b* 1 /2 . 53 ) ; 
printf ("\n%s  %d"," Recommended  number  *  " ,$T); 
printf  ("\n%s", "Want  to  change  ?  (Y/N) " ) ; whi 1 e( I kbhi t () ) ; 
res=getch () ; res=tol ower (res) ; i f  (res== ' y ' ) 

{printf  ("\n%s" , "Input  number  of  spatial  nodes  "); 
scanf  ("%d",&ST);} 

settextposition  (16,30);  settextcolor  (28);  outtext( "Working 

)  j  - 

re= (q*904000/ (d*tg) ) *pow(298/tg , .75); 

cpg«=(6.449+.001413*tg-(8.07e-08)*pow(tgl2))/29; 

pr=cpg/(.0855+cpg) ;nu=.023*pow(re, .793)*pow(pr, .333) ; 

ka=.0001845*cpg*pow(tg/298, .75)/pr; 

hfc=ka*nu/d; 

delx“(float)  1/ST; 
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iter 


de1t=(f1oat)  pow(delx,2)*rho*cp*pow(l ,2)/(ks*36G0) ; 
nz*(int) (1+n/delt) ; 
if  (3FL[l]==ll){for  (j=0; j<-nz+4; j++) 

■'for  (i=0; i<=ST; i++)  {th[ij  [j]=ti/tg; }}} 
if  (BFL[2]--=11)  {for  (j=0;j<=nz+4;j++){for  (i=0;i<=ST;i++){ 
x=(float) (ST-i )/ST;th[i] [j]=(x*ti/tg)+(l-x)*ta/tg;}}} 
if  (BFL[4]==11)  {  ts=(ks*ti*lp+kp*ta*l)/(kp*l+ks*lp) ; 
for  (j=0; j<=nz+4; j++){for  (i=0;i<=ST;i++){ 
x=(float) (ST-i ) /ST; th Li ] [j]=(x*ti/tg)+(l-x)*ts/tg;}}} 
if  (BFL[3]==ll)  {for  (j  =0; j<=nz+4; j++) {for  (i=0;i<"ST;i++){ 
x=(float)  (ST-i)/ST;th[i][j]=(x*ti/tg)+(l-x)*ts/tg;}}} 
switch(cres){  /*units  of  cs  and  eg  not  in  right  dimension*/ 
case  'a1: 

for  (j=0;j<=nz+4;j++){for  (i=0;i<=tST;i++)  { 
cs[i][j]=coav;}} 
break; 
case  *  b ' : 

for  (j=0;j<=nz+4; j++) [for  (i=Q;i<=ST;i++){ 
x=(float) (ST-i )/ST;cs[i] [j]=coav*x*2;}} 
break; 
case  1 c ' : 

for  (j=0;j<=nz+4;j++){for  (i=0; i<=ST; i++){ 
x= (float) (ST  i)/ST;cs[i] [j]=2*coav*(l-x) ;}} 
break; 

j j=0; k=12; 

:  dth=th  [0]  [j  j+1] ;  for  (j=jj; j<=k+2;j++){ 

hr=(5.01e-13)*pow(tg,3)*(l-pow(th[0] [j+l] , 4) ) / (1-th [0] [j+1]); 
gr=fabs(((2.09el4)*pow(d,3)/pow(tg,3.5))* 

(1-th [0] [j+1])/ (l+pow(th[0j [j+1] ,4.5))); 
hc=.525*(ka/d)*pow(pr*gr, .25) ; 
ht=sqrt(pow(hfc, 2)+pow(hc,2))+hr; 

“  '  [j+l]=.5*(th[l]  [j]+th[l]  [j+l]+(ht*l /ks) * 

[0]  [j]-th[0]  [j+lj  )*delx) ; 

'  "=-ll){  th  [ST]  [j+1]  =th  [0]  [j+1] ; } 

==ll){gra=fabs((2.094el4)*(pow(d,3)/pow(ta,3.5))* 
j+1])/ (l+pow(th[ST] [j+1] ,4.5))); 
ha=.525*(ka/d)*pow(pr*gra, .25) ; 


th  [0] 
(2-th 
if  (BFL[1 
if  (BFL[2 
(l-th[ST] 


th [ST] [j+l]=.5*(th [ST-1] [j 
(ha*l/ks)*((2*ta/tg)-th[ST' 
if  (BFL[3]==ll){th[ST] [j+1] 

(kd*l/ks)*pow(((j+j+l)*.5*delt*ad) ( 
(2*(ts/tg)-th[ST] [j]-th[ST] " 
if 


+th  [ST-1]  [j+1]  + 

Cj] -th  [ST] [j+l])*delx) ;} 
=.5*(th[ST-l] [j]+th [ST-1] [j+l]+ 
.5)* 

[j+1] )  *del  x) ;} 

=.5*(th[ST-l] [j+l]+th[ST-l] [j]+ 
][j+l])*delx);} 


(BFL[4J==ll){th[ST] [j+1]  .  . 

(kp*l/(lp*ks))*(2*(ta/tg)-th[ST][j]-th[Sl 
for  (i*l;i<*ST-l;i++){ 
th[i][j+l]  =  .25*(th[i+l][j+l]+th[i+l]  [j]+th  [i-1]  [j]+th  [i -1]  [j+1] ) ; }} 
tt=dth-th [0] [j j+1] ; 
if  (fabs(tt)>=le-07){  goto  iter;} 
k=k+l ; j j=j j+1 ; i f  (jj<=nz+4){goto  iter;} 
output () ; 

if  (resp[0]=='m' ) {goto  start;} 
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output () 

{ 


(BFL[7]«-11) 
(BFL[6j— 11) 


if  (BFL[7 
if  (BFL[6_ 
return (0); 


11)  store(); 
11)  plot(); 


menu  () 

input:  _clearscreen  (  GCLEARSCREEN)  ,-^set videomode  (_DEFAULTMODE) ; 

_settextcolor  (147;_settextpositTon  (1,30); 
outtext  ("Hot  Gas  Properties"); 

<f*60 .96*he*wi / (he+wi ) ;_$ettextcol or (2) ; 
_settextposition(2,5);_outtext  ("Gas  (F)low  Rate  »  "); 
sprintf  (bufl,"%4.0f",q);_outtext(bufl);  outtext  ("  cfm"); 
_settextposition(2,50) ;  outtext  ("(G)as  Temperature  *  "); 
sprintf  (bufl, "%3.1f" , (tg-273)*l .8+32) ;_outtext(bufl) ; 

_outtext  ("  F\n");_settextcolor(14);_$ettextposition  (3,30); 

_outtext  ("Wall  Properties");  settextcolor (2);  settextposition(4,5); 
_outtext  ("(C)onductivity  ■  "7;sprintf(bufl,"%5.2f",ks*241) ; 
_outtext(bufl);_outtext("  Btu/ft-hr -F") ;  settextposition(4, 50) ; 
“outtext  ("(S) pacific  Heat  «  ");sprintf(Fufl,"%1.2f",cp); 
_outtext(bufl) ;_outtext("  Btu/lb-F") ;  settextposition(5,5) ; 

“outtext  ("(D)ensity  *  ");  sprintf (bufT,"%3. Of ",(rho/0. 016)); 
_outtext(bufl);  outtext("  1b/ft^3");  settextposition (5, 50); 

“outtext  (“(T)hTckness  *  ");sprintf(Eufl,"%2. If", (1/2.54)); 
_outtext(bufl);  outtcxt("  in.");  settextposition(6,5) ; 

_outtext  ("(H) eight,  Width  *  ");“ 

sprintf (bufl,"%2. If  %s  %2.1f",he,"x",wi);_outtext(bufl); 

_outtext("  ft.");  settextposition(6,50) ; 

“outtext  ("(I)nitTal  Temperature  a  "); 

spri ntf (buf 1 , "%3 . If" f ( (ti -273) *1 .8+32) ) ;_outtext (buf 1) ; 

_outtext("  F");  settextcolor(l4) ;_settextposition  (8,30); 

“outtext  ("Boundary  Conditions") ;“settextposition  (9 . 10) ; 
“settextcolor  (BFL[1)) ;_outtext  (¥(1).  Both  sides  heated"); 
“settextposition  (9,50)7_settextcolor  (BFL[2]); 

_outtext  ("(2).  Air  cooled  side");  settextposition  (10,10); 
^settextcolor  (BFL[3] ) ;  outtext  ("13).  M^st  soil  side"); 
“settextposition  (10,507;_settextcolor  (BFL[4]); 

“outtext  ("(4).  Insulated“side") ;_settextposition  (12,30); 
Jsettextcolor  (14);  outtext  ("Output  Choices"): 

'settextposition  (if,  10);  settextcolor  (BFL[5j); 

If  (BFL[51bb3){  outtext  (*(5).  IBM  Printer");} 
if  (BFL[5J**ll)X_outtext  ("(5).  Epson  Printer');} 
_settextposition”(13,30);_settextcolor  (BFL[6]); 

“outtext  ("(6).  Plotter")T_settextposition  (13,50); 

“settextcolor  (BFL[7]);  outtext  ("(7).  Disk"); 

“settextposition  (14,307;  settextcolor  (14); 

“outtext  ("Chemical  Agent^) ;_settextposition  (15,10); 
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settextcolor  (BFL [8] > t  outtext  ("(8).  HD");  settextposition  (15,30); 
_settextcolor  (BFL [9j  "outtext  ("(9).  GB")  .‘“settextposition  (15,50); 
settextcolor  (BFLLO];,  outtext  ("(o).  VX"); 

Tf  ( ( BK. [8]  +BFL [9]  +BFL [0] ) ==1 7 ) { 
if(  BFL[6]+BFL[73)>=14){ 
if((BFL[l]+BFLh3+BFL[3t+BFL[4])==20){yp=5; 
settextposition  (16,30);_settextcolor  (12); 
outtext  ("(R)eady  to  Calculate");}}} 
printf  ("\n\n%s  V'Set  up 

problem  by  pressing  letter  or  number  key  for  each  "); 
printf  ("\n%s\n", "variable  you  wish  to  change  or  set"); 
wh i 1 e ( ! kbh i t () ) ;  resp [0] =getch ( ) ; resp [0] =tol ower ( resp [0] ) ; 
switch  (resp[0]){ 
case  ' g ' : 

printf  ("  Input  gas  temperature  (F)  "); 
scanf  ("%f",&tg); 
tg=((tg-32)/1.8)+273; 
break; 
case  ' f' : 

printf  ("Input  volumetric  flow  rate  (cfm)  "); 
scanf  ("%f",&q); 
break; 
case  ' s' : 

printf  ("Input  specific  heat  (Btu/lb-F)  "); 
scanf  ("%f",&cp); 
break; 
case  ' i ' : 

printf  ("fnput  initial  temperature  (T)"); 
scanf  ,  .  ,&ti) ;tia(.5555*(ti-32)+273); 

break; 
case  ' c ' : 

printf  ("Input  Thermal  conductivity  (Btu/hr-ft-F)  "); 
scanf  ("%f",&ks);ks-ks/241; 
break; 
case  ‘d‘ : 

printf  ("Input  density  ( 1 b/f t^3)  ”); 
scanf  ("%f",&rho) ;rho=rho*.016; 
break; 
case  ' t 1 : 

printf  ("Input  wall  thickness  (in.)  "); 
scanf  ( "%f",&l); 1=1 *2.54; 
break; 
case  ' h ' : 

printf  ("Input  room  height  (ft)"); 
scanf  ("%f",&he); 
printf  ("Input  room  width  (ft)"); 
scanf  ("%f",&wi); 
d=30,48*(4*he*wi/(2*he+2*wi)); 
break; 
case  ' r' : 

if  (yp<=l){_settextcolor(12) ;_settextposition(20, 1) ; 
_outtext("You  have  not  selected  the  minimum  blue  items 
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case 


case 


case 


case 


—press  enter  to  continue"); 
while(lkbhit());break;} 
r^  *urn  (0) ; 

1  1  1  ° 

for* (ia2;i<*4:i++){BFL[i]a3;} 

BFL[13»14-BFL[1]; 

break; 

l  2'  ; 

BFL[2]a14-BFL[2] ;BFL[l]a3; 

*'FL[3]-3;BFL[43-3; 
intf  ("Input  air  temperature  (F)  ) ; 
anf  ("%f" ,&ta);ta»273+.5555*(ta-32) ; 
iireak; 

'4' : 

BFL[4l“14-BFL[4] ; 

BFL 1 1 J “3 ; BFL [2] =3 ; BFL [3] *3 ; 

printf  ("%s\n\r", "Thermal  conductivity 

of  insulation  -  0.05  Btu/hr-ft-F  ;  Want  to  Change  (Y/N)"); 

whi 1 e( I kbhi t ( ) ) ; resagetch () ; res«tol ower (res) ; 

if(res**'y')tprintf("Input  conductivity  of 

insulation  (Btu/hr-ft-F)"); 

scanf  ("%f" ,&kp) ;kpakp/241;} 

printf  ("Input,  thickness  of  insulation  (in)"); 

scanf  ( "%f " , &1 p) ; 1 p*l p*2 .  54 ; 

printf  ("\n\r%s\ "Input  outside  air  temperature  (F)'); 

scanf  ("%f" ,&ta) ;ta*273+.5555*(ta-32) ; 

break; 

•  3  • . 

BFL[3]=14-BFL[3];8F1.[4>3; 

for  (i=l;i<=2;i++){(BFL[ij»3);}  . 

printf ("Soil  temperature  a  55  F;  Want  to  change  ?  (Y/N)  ) 
while  (IkbhitQ)  ;resagetch()  ;res=to  lower  (res)  ; 
if(reS**'y'){printf  (“\n\r%s"," Input  soil  temperature  (F)" 
scanf  ("%f" ,&ts) ;ts°273+.5555*(ts-32) ;} 
printf  ("\n\r%s  %2.1f  %s", "Moisture  content  a  ",wa,"% 

Want  to  change  7  (Y/N)  "); 

while  (!kbhit());res  agetch();resa  tolower(res) ; 
if  (res»*y)  {printf  ("\n\r%s" , "Input  soil  moisture 
content  (%) " ) ; 

$canf  |  '•  j  •  1 

rhod* (wa/ 100) + ( 1 -wa/ 100) * ( 1 . 5) ; cpda (wa/ 100) + . 2* ( 1 -wa/ 100) ; 
kd*wa*.000014+.0024*(l-wa/100) ;adakd*3600/ (rhod*cpd);ta»ts 


) 


break; 
case  '7' : 

Bfl[7] =14-BFL[7] ; 

printf  ("(S)torenew  data  or  (R)etrieve  old  data  ?"); 

while(lkbhit());res=getch();res=tolower(res); 

printf  ("\n\r%s\" Input  file  name  "); 

scanf  ("%s" ,fname) ; 

if  (resaa'r'){retrieve(); 

BFL[7]a14-BFL[7] ;} 
break; 
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case 

case 

case 

case 

case 


'5* : 

BpL[53*i4>  HFL[5]; 
break; 

'6': 

BFL[6]a14-BFL[6] ; 
break; 

•8': 

BFL [81 - 1 1 ; BFL [9] «3 ; BFL [0] -3 ; 

chem(); 

break; 

•  g  •  . 

BFL [9] * 1 1 ; BFL [8] -3 ; BFL [0] =3 ; 

chem(); 

break; 

•O': 


BFL  COl *1 1 ; BFL [8] =3 ; BFL [9] =3 ; 
chem(); 
break; 
default: 

printf  ("Try  again");  } 
goto  input; 


} 

chem() 

clearscreen  (_GCLEARSCREEN);if(  setvi^omodefERESCOLOR^Q) 
Xypix=200;if(_setvideomode  (  HRE516C0L  }*=0){  setvideomode(_KRESBW) ;}} 
setcolor(ll);  rectangle(  GB0RDER,2O,O,i7O, (.2557*ypix)); 
rectangl e (_6B5RDER , 210 , 0 , 370 , ( . 2857*ypi x) ) ; 

"rectang 1 e (  GBOROER ,410,0.570, (. 2857*ypi x) ) ; 

_setcolor (1?) ;  jnoveto (20 , 50) ;  1 i neto( 170 , 50) ; 
jnoveto(210,0)7  lineto(370, (.2857*ypix)); 

_moveto(410, (.2557*ypix)) ;_lineto(570,0); 
_settextcolor(12);_settextpositi m(4, l);_outtext("C") ; 

^settextpos i ti on (9712) ;_outtext ( "x" ) ;_settextpos i t i on (9 , 37 ) ; 
js  uttext ( " x " ) ;  _set text pos i t i on ( 9 , 62 ) ; "out text ( " x " ) ;  set textco 1  or ( 10 ) ; 
7$ettextposi ti on (2,12) ;_outtext ( " (a) "7 ;_settextposi ti on (2,37); 

^outtext ("(b)");  settextpos i t i on (2 , 62 ) ;~ou ttext ("(c)") ; 
_settextcolor(llj;_settextposition(12, lj; 
j>uttext(" Select  a~b  or  c  as  initial  agent  distributions"); 
wh i 1 e ( I kbhi t ( ) ) ; cres*getch () ; cres=tol ower (cres) ; 
settextpos i t i on ( 14 , 1 )  ;_outtext ( 

TInput  average  aaent  loading  (mg/g)  "); 
scanf ( "%f " , &coav) ;  settextposi tlon (15,1); 

_outtext(" Select  (E)vaporation  or  (P)yrolysis  as  release  mode  "); 
wh i 1 e ( I kbh i t ( ) ) ; mode=getch ( ) ; mode*tol ower (mode) ; 
return (0) ; 
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int  z; 

clearscreen  (J3CLEARSCREEN) ; if (  set v i deomode ( JERESCOLCR) ==0 ) 
7ypixa200; if (_setvi deomode  (^HRESlBCOLCRj^Q) 

{_setvi deomode (  HRESBW);}}  setcolor(ll) ;  setcolor(ll) ; 
rnoveto  (0,3);  Tineto  (0,07;  jnoveto  (0,3j;_lineto  (7,3); 

Suffer  *(char  Tar  *)malloc((unsigned  int) 

_i mages ize{0, 0,1,1)); 

_getimage(0,0, 10, 10, buffer) ; 

replot :  clearscreen  (_GCLEARSCREEN) ; i f (_setvi deomode (_ERESC0LOR) ==0) 
{ypTx*200;if (_setvi deomode  (_HRES16tOLOR)=sO) 

_setvi deomode (_HRESBW) ; }}_setcol or (1 1) ;_settextcol or ( 1 1) ; 
printf ("Select  variables  to  be  plotted") ;_settextposition (2, 1) ; 
outtext  ("Plot  (a)  Temperature  (b)  Cone,  of  adsorbed  phase 
Tc)  Gas  phase  cone.  "); 

_settextposition(3,l) ;_outtext("  (d)  Fraction  of  agent  left  ?"); 
while(lkbhitO)  ;cres=getch()  ;cres*tolower(cres); 
if(cres=*'d' ) {presp [OJ^'t1 ;} 
if(cres!*'d,){  settextposition(4,l); 

_outtext  ("(T)Ime  or  (D)epth  Profile  ?"); 

while  (I kbhi t ( ) ) ; presp [0] =getch () ; presp [0] =tol ower( presp [0] ) ; } 

switch  (presp [0J){ 

case  ' t ' : 

if  (crest'd1 )  {break; } 

settextposition(4,l);  outtext("Available  Depths  ="); 

Tor  (j=0; j<=ST;j++){prIntf  ("%4.2f  " ,j*l*delx/2.54) ;} 
rcoord*jjettextposition();rcoord.row  += (ST/ 12) 

_settext pos i t i on ( rc  oord . row , 1 ) ; 

_outtext(" Number  of  Depths  to  plot  *  ?  "); 

SCanf  { 11  %d &NP)  * 

for  (i=l;i<=NP;i++) {printf  ("%s  %d  %s","Depth  #",i,"  *  "); 
scanf  ("%f" ,&pl [i]) ;} 
break; 
case  ' d  * : 

settextposition(4, 1) ;  outtext("Available  Times  ="); 

Tor  (j=0;j<=nz;j++) {printf  ("%4.2f  ",j*delt);} 
rcoord=_gettextpos i ti on (); record . row  += ( nz/ 12) +2 ; 

_set textpos i t i on ( rcoo  rd . row , 1 ) ; 

_outtext(" Number  of  Times  to  plot  =  ?  "); 
scanf  ("%d",&NP); 

for  (i=l; i<=NP; i++) {printf  ("%s  %d  %s","Time  #",i,"  =  "); 

scanf  ("%f " , &pl [i] ) ; > 

break; 

_clearscreen  (_GCLEARSCREEN) ;  setcolor(ll) ; 

_rectangle  (  GBQRDER, 110, 1,63$, .8*ypix);ovr();  settextcolor(ll) ; 
_settextposItion(24,l) ;_outtext(" (p)lotter  (rjeplot  prn(s)crn 
prin(t)  (m)enu  or  (e)xit  ?"); 

_settextcolor(12) ; 
switch  (cres){ 
case  'a': 

tt=(tg-273) *1.8+32; 

for  (yp=0;yp<=.95*tt;yp+att/10){ 


settextposition  (20.3* ( 1-yp/tt) ,9) ;sprintf (bufl, "%4.0f" ,yp) ; 
“outtext (bufl) ;}_settextcolor(14) ;_settextposition  (10,1) ; 
“outtext  ("Temp,“F"); 

Freak; 
case  1 b ' : 

tt=4*coav; 

for(yp=0;yp<=.95*tt;yp+att/10){ 
settextposi ti on (20 .3* (1-yp/tt) , 10) ; spri rtf (bufl , "%2 . If " ,yp) ; 
“outtext (bufl ) ; }  settextcol or ( 14) ; _settextpos i t i on ( 10 , 1 ) ; 
outtext ( "C(ad) fT) ; 

"settextposi tion (11, 1) ;_outtext("mg/g") ; 
break; 
case  1 c 1 : 

tt=4*coav; 

for  (yp=0;yp<s!.95*tt;yp+stt/10){ 
settextposition  (20.3*(l-yp/tt) , 10) ; 
spri ntf (buf 1 , "%2 . If” ,yp) ;_outtext (bufl) ; } 

_settextcol or (14) ;  settextpos i t i on ( 10 , 1 ) ;_outtext ( "C (g) t " ) ; 

settextposition(lI, 1) ;_outtext("mg/L") ; 

Freak; 
case  ' d ' : 

for  (yp=0;yp<=.95;yp+=.l){ 

_settextposition(20.3*(l-yp) , 10) ; sprintf (bufl , "%1 .2f " ,yp) ; 
_outtext(bufl) ;} 

_settextcolor(14) ;  settextposition(10, 1) ;  outtext("Fraction"); 
_settextposition(lI,  1)  ;_outtext("of  agent17)  ; 

settextpos i t i on ( 12 , 1 ) ; _outtext ( " 1  eft " ) ; 

Freak; 

switch  (presp[0]){ 
case  ’ t 1 : 


dths:nz*delt  ;_settextcolor(12) ; 

for  (tt=0;tt<=l;tt+=.l){  settextposition  (22,14+65*tt); 
spri ntf (buf 1 , "%2 . If " , tt*3th) ;^outtext (buf 1 ) ; } 
_settextcolor(14) ;_settextposTtion  (23,40) ; 

_outtext  ("Time,  hours") ;_setcolor(14); 
switch (cres){ 
case  ' d ’ : 

moveto( 110,0) ; 

Tor(j=l;j<=nz;j++){yp=fnpl(); 
tt=110+(float)529*j/nz;_1ineto(tt,yp) ;} 
break; 
default: 

settextposition(2,75) ;  outtext("Depth") ; 

?or  (jj=l;jj<=NP;jj-H-)-'“ 


j.l+pl [jj ]7 (l*delx/2.54) ) ; j=0; 
moveto  (i: 


i=(int) 

!ll0,fnpl ());_setcolor(15-jj) ; 
’settextcolor(15-jj) jsprintf 
Tbufl , "%2.2f",i*(l*delx/2.54)); 
settextposition  (jj+2,75);  outtext (bufl) ; 
Tor  (jsl;j<=nz;j++) (yp=fnplT) ;  % 

tt=110+(float)529*j/nz;_lineto  (tt,yp) ;}} 
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break; 

^or  (i  s  163;  i<=640;  i  +B  53) { 
putimage  (i , (.8*ypix-3) , buffer,  GOR);} 
tta(tg-273) *1.8+32; 

for  (yp=0;yp<= . 95*tt ;yp+=tt/ 10) {?* ( .8*ypi x-3) * ( 1-yp/tt) ; 
putimage  (110,z,buffer,_GOR) ;} 

Freak; 
case  *  d 1 « 

dth°l/2.54;  settextcolor(12); 

for  (tt=0;tt<=.99;tt+=(.l)){_settextposition  (22, 14+65*tt) ; 
sprintf (buf 1 , "%2.2f " , tt*dth) ;_outtext(bufl) ;} 
_settextcolor(14) ;_settextposition  (23,37); 

_outtext  ("Distance,  inches") ;_setcol or (14); 

settextposition(2,75) ;  outtext ("Times") ; 

Tor  (jj*l;jj<“NP; jj++){J=(int)  (.1+pl [jj]/delt) ;i*0; 
moveto  (110 , fnpl () ) ;  _setcolor(15-jj);  settextcolor(15-jj) ; 
spri ntf (buf 1 , "%2 . 2f " , j *del t ) ; 
settextposition(jj+2,75) ;  outtext (buf 1) ; 

Tor  (i=l;i<='ST;i++){yp=fnpT(); 

tt=110  +(float)529*i/$T;_lineto  (tt ,yp) ; }} 

for  (i  *  163;  i<=640;  i  +=  53) { 

_putimage  (i , (.8*ypix-3) , buffer, _G0R) ;} 

tt*(tg-273) *1.8+32; 

for  (yp=0;yp<=.95*tt;yp+=tt/10){ 

z»(.8*ypix-3)*(l-yp/tt);_putimage  (110, z, buffer, _G0R) ;} 
break; 

} 

out:  while(!kbhit());resp[0]agetch();resp[0]a,tolowe»  (resp[0]); 
switch  (resp[0J){ 
case  's': 

prnscrn(); 
goto  out; 
break; 
case  ' p* : 
hp(); 
goto  out; 
break; 
case  ' t 1 : 

print(); 
goto  out; 
break; 
case  ' r* : 

_setvideomode(J)EFAULTMODE) ; 
goto  replot; 
break; 
case  'e': 

goto  ret; 
break; 
case  'ml : 

goto  ret; 
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break; 

default: 

printf  ("program  complete");  } 
goto  replot; 

ret:  jsetvideomode  (_DEFAULTMODE) ; 
return  (0) ; 


switch  (cres){ 
case  'a': 

yp  »  .8*ypix*(tg*(l-th[i] [j])/ (tg-255) ) ; 
break; 
case  ' b ' : 

yp  *  .8*ypix*(l-cs[i][j]/(4*coav)); 
break; 
case  ’ c  * : 

yp  *  .8*ypix*(l-cg[i][j]/(4*coav)); 
break; 
case  1 d 1 : 
yp=0; 

for(ia0;i<*ST-l;i++){ 

yp*yp+.5*(cs[i][j]+cs[i+l][j]+cg[i][j]+cg[i+l][j])/ST;} 

yp=.8*ypix*(l-yp/coav) ; 

break; 

return (”%f",yp); 

jp() 

stream  «  fopen  ("COMl'Y'wr"); 

setvbuf (stream, buf,  lOFBF.sizeof (buf) ) ; 

f printf  (stream,  "\33%s",".N;19:M  ); 

f printf  (stream, "\33%s" , " .H180; ;17:") ; 

f printf  (stream," IN; SP0;VS25;"); 

fprintf  (stream, "IP1700, 1520, 8175, 7350; SP1;"); 

f printf  (stream,  "PUPA1700, 1520; PDPA1700 , 1520, 1700,7350,8175, 

7350,8175, 1520,1700,1520;") ; 

fprintf  (stream,  "PDPA1720, 1520, 1720,7350,8175,7350,8175, 1540, 
1720, 1540; PU;"); 

fprintf  (stream,  "Si . 15 , . 25; tl 2 ,0; sp2; " ) ; hand() ; 

switch (cres){ 
case  'a' : 

fprintf  (stream, "%s\3\  "SPl;DI0,l;PUPA1700,3737;CP-7,3; 

LBTemperature,  F");hand(); 

for  (d“0;d<ai1.05;d+“.l){ 

tt= (255. 3+d* (tg-255. 3) )/tg; 

hy=  1520+4309* (tt-255 .3/tg) /( 1-255 ,3/tg) ) ; 

hx»(32+1.8*(tg*tt-273)); 

fprintf  (stream, "%s  %d  %4.0f  ","TL2,0;DI;PUPA",1700,hy); 
fprintf  (stream, "%s  %4.0f\3",";YT;CP-5,0;LB",hx);hand(); 
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if (d<l){hy=hy+216;fprintf (stream, "%s  %d  %4.0f  %s" ,  "TLl.O; 

PUPA" ,1/00, hy, “ ; YT") ; }} 

break; 

C3Se  fprintf  (stream, "%s\3",  "SP1 ; DIO , 1 ; PUPA1700 , 3737; CP-14,3 j 
LBConcentration  in  Solid,  mg/g");hand(); 
for  (d“0;d<=*1.05;d+st.l){ 
hy»(1520+4309*d); 
hx=d*4*coav; 

fprintf  (stream, "%s  %d  %4.0f  " ,"DI; PUPA" ,1700, hv); 
fprintf  (stream, "%s  %3.1f\3\";YT;CP-5,0;LB",hx);hand();} 
break; 

CaSe  fprintf  (stream, "%s\3\  "SP1;DI0, 1; PUPA1700, 3737; CP-15; 
LBGas  Phase  Concentration,  mg/L") ;hand() ; 
for  (d*0;d<=!I.05;d+*.l){ 
hy»(1520+4309*d); 
hx*d*4*coav * 

fprintf  (stream, M%s  %d  %4.0f  “ , " DI ; PUPA" , 1 700 . hy ) ; 
fprintf  (stream, ”%s  %3. lf\3“ , “ ;¥T;CP-5,0;LB" ,hx) ; handO ; > 
break; 

fprintf  (stream, "%s\3" ,  "$Pl;DI0,l;PUPA170O,3737;CP-ll,3; 
LBFraction  of  Agent  Left");hand() ; 
for  (dst0;d<*1.05;d+*.l){ 

fprintf+istream/%s'%d  %4.0f  »,"TL2,0;DI; PUPA", 1700, hy) ; 
fprintf  (stream, "%s  %1. 2f\3 \\-YT;  CP-5,0;  LB"  ,hx) ;  hand  () ; 
if (d<l){hy=hy+216;fprintf (stream, "%s  %d  %4.0f  %s", "TL1.0; 
PUPA", 1700 ,hy,";YT");}} 

i  i 


} 


break; 


switch  (presp[0]){ 
ca$e  • ^  ♦ 

for*  (d=0;d<='(float)1.05*n;d+«l+(int)n/10){tt=d; 

hx» ( 1700+6475* (fl oat) d/n); 

fprintf  (stream, "%s  %4.0f  %d  %s  %2.1f\3  , 

nTL2,0;DI;PUPA" ,hx, 1520," ; XT; CP-1, -1; LB" ,tt); hand () ;} 

for  (d*0;d<*(float)n;d+*.1667*(l+(int)n/10)){ 
hx«1700+6475*(float)d/n; 
f printf (stream, "%s  %4.0f  %d  %s", 

" TL1 , 0 ; PUPA" , hx , 1 520 , " ; XT ; " ) ; hand ( ) ; }  icon  ro  . 

fprintf  (stream, "%s\3" , "SP1; SI; PUPA4937, 1520; CP-6, -3; 
LBTime,  hours") ; hand () ; 
switch (cres){ 

case  1  d1:  . 

f pr i ntf ( stream , " S  P2 ; PUPA1 700 , 5829 ; " ) ; 
for(j=i; j<»nz; i++) (hy=fnhp() ; 
hx=»1700+  (f  1  oat)  647  5*j  /nz; 

fprintf (stream, "%s  %4.0f  %4.0f","PDPA",hx,hy);hand();} 
break; 
default: 
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for  (k=l;k<=NP;k++){i=(int) (.1+pl [k] / ( 1 *del x/2 . 54) ) ; 
j=0;hy*fnhp(); 

fprintf  (stream, "%s  %4.0f","SP2;PUPA1700",hy); 
for  (d=0 ; d<=n/del t ; d++) { j  « ( i nt ) d ; 
hx=(1700+6475*(j*delt)/n) ;hy=fnhp() ; 
fprintf  (stream, "%s  %4.0f  %4.0f", 

" PDPA" , hx , hy ) ; hand () ; } 
fprintf  (stream,"%s  %2.1f  %s\3", 

"PU;DI;SI.10, .15;CPl,0;LB",pl [k] ,"in.");} 

ilreak; 
case  1 d ' : 

for  (d=0; d<=l .05; d+= (2.54/1)) {tt=d*l /2 . 54; hx= ( 1700+6475*d) ; 
fprintf  (stream, "%s  %4.0f  %d  %s  %2.1f\3", 

"TL2  f  0;  DI ;  PUPA",  hx,  1520, 11 ;  XT;  CP-1, -1;  LB",  tt);  hand  ();} 
for  (d=0;d<*l;d+-.'.54/ (2*1 ) ) {hx=1700+6475*d; 
fprintf (stream, "%s  %4.0f  %d  %s" , 

"TL1 ,0; PUPA" , hx , 1520 , " ; XT; “ ) ; hand () ; } 
f pri ntf  ( stream , "%s\3 " , "SP1;SI; PUPA4937 , 1520; CP-6 , -3 ; 

LBDepth,  inches"); hand (); 

for  (k-1 ; k<=NP; k++) { j  =  ( int) ( . 1+pl [k] /del t) ; i =0; hy^fnhp () ; 
fprintf  (stream, "%s  %4.0f  %s  %2.1f  %s\3“, 

"SP2;PUPA1700",hy, ";DI;SI.10, .15;CP.5,-.5;LB< — ",pl [k] , "hr") ; 
hand() ; 

fprintf  (stream, "%s  %4.0f","PUPA1700",hy); 

for  ( i =0 ; i<=ST ; i ++) {d= ( f 1  oat) i /ST ; hx= ( 1 700+6475*d) ; hy =f nhp ( ) ; 

fprintf  (stream, "%s  %4.0f  %4.0f","PDPA",hx,hy);hand();}} 

fprintf  (stream, "%s\3","SPl;PUPA4937,7350;Di ;Si;CP-13,-l; 

LBHotgas  Desorption  Process"); 

fprintf  (stream, "%s  %d\3" ,"SP2;PUPA1700,7350;CP2,-2;Si .15, .25; 

LB  Number  of  spatial  nodes  *  ",ST);hand(); 
hy=(tg-273) *1.8+32; fprintf  (stream, "%s  %4.0f  %s\3","CP; 

LB  Gas  Temperature  =  ",hy,"F"); 

fprintf  (stream, "%s  %5.0f  %s\3",“CP;LB  Gas  Flow  Rate  *  ",q,"cfm"); 
fprintf  (stream, "%s  %2.1f  %s\3","CP;LB  Wall  Thickness  =* 

",  1/2. 54, "in."); 

fprintf  (stream, "%s  %2.0f  %s  %2.0f  %s\3","CP; 

LB  Dimensions  =  ",he,"  x  " ,wi , ‘ ft.") ; 
if  (BFL[1]*=11) {fprintf  (stream, "%s\3" , "CP; 

LB  Both  sides  of  wall  heated");} 
hx=  ( 1 .8* (ta-273)  +32) ;  i f  (BFL[2]«11){ 
fprintf (stream, "%s  %4.0f  %s\3","CP; 

LB  Outer  wall  air  cooled  with  air  temperature  =  ",hx,"  F");} 
hx=l .8*(ts-273)+32; i f  (BFL[3]=»11){ 
fprintf  (stream, "%s\3", "CP; 

LB  Outer  wall  in  contact  with  moist  soil"); 
fprintf  (stream,  "%s  %2.0f  %s  %4.0f  %s\3","CP; 

LB  Water  content  -  ",wa,"%  T  =  ",hx,"F");} 
if  (BFL [4] ==1 1 ) {fpri ntf  (stream, "%s  %2.0f  %s\3","CP; 

LB  Outer  wall  insulated  with  " , lp/2.54, "in.  @k  ■  0.05  Btu/hr-ft-F");} 
fprintf  (stream,  "PUPA0,0;SP0") ; 
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} 


fflush  (stream); 
fclose; 
return (0) ; 


store() 
{ 


} 


stream  =  f open  (f name,  V'); 

fprintf (stream, "%d  %f  %f  %f  %f  %d  *f  %f 

,ST,tg,ti , q, 1 ,he,wi ,nz, lp,ta,delx,delt,coav,cres); 
for  (i=0;i<=9;i++) {fprintf  (streamtl'%cl  \BFL[i]);} 

for  (1=0; i<sST; i++|  'fprintf (stream, "%f  \th[i][j]);}} 
fclose  (stream); 
return (0) ; 


%c  " 


retrieve() 

{ 


lieader() 


stream  *  fopen  (fname,  "r");  ...  .. 

fscanf(stream,"%d  %f  %f  %f  %f  %f  %d 

,  &ST ,  &tg ,  &t  i ,  &q ,  &1 ,  &he ,  &wi ,  &nz ,  &1  p ,  &ta ,  &del  x ,  &del  t ,  &coav ,  &cres) ; 

for  (i=0;i<*9;i++){fscanf  (stream, “%d  “ ,&BFL[i]) ;} 

lor  (isa0;i<!,ST;i++|  fscanf (stream, "%f  “ ,&th[i] [j]) ;}} 
fclose  (stream) ;BFL[7]=S3;BFL[6]=11; 
output () ; 
return (0) ; 


He 


fprintf  (stream,  "  Hotgas  Desorption  Process"); 
switch (cres){ 
case  'a1: 

fprintf  (stream, "\n\r%s",  "Temperature,  (F)  ); 
break; 
case  1 b ' : 

fprintf  (stream,  "\n\r%s\ 

"Adsorbed  Phase  Concentration,  (mg/g)  ); 


break; 

CaSe  fprintf  (stream, "\n\r%s\  "Gas  Phase  Concentration,  (mg/L)"); 
break; 

fprintf  (stream, "\n\r%s  %d",  "  Number  of  spatial  nodes  =  \ST); 
fprintf  (stream, "\n\r%s  %4.1f  %s", 

"Gas  Temperature  *  " ,  ((tg-273)*1.8)+32, '  F  ) ;  „  f  ,,v. 

fprintf  (stream, "\n\r%s  %4.1f  %s",  "Gas  Flow  Rate  -  ,q,  cfm  ), 

fprintf  (stream, "\n\r%s  %2-lf  %s",  "Wall  Thickness  =  ,1/2.54,  in.  ), 

fprintf  (stream, "\n\r%s  %2.1f  %s  %2.1f  %s  , 

"Dimensions  =  ",he,"  x  ",wi,"ft."); 

if (  BFL[1]*-11)  {  fprintf  (stream, "\n\r%s" , 

"Both  sides  of  wall  heated");} 
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if(  BFL[2]*=11)  {  fprintf  (stream, "\n\r%s  %4. If  %s", 

"Outer  wal 1  air  cooled  with  air  temperature  =  ", (ta-273)*1.8+32, "  F");} 
1f(  BFL[3]“11)  {  fprintf  (stream, ,t\n\r%s" , 

"Outer  wall  in  contact  with  moist  soil");} 

if(  BFL[4]=*11)  fprintf  (stream, "\n\r%s  %2.1f  %s  %1.3f  %s", 

"Outer  wall  insulated  with  ",lp/2.54,"in.  of  insulation 
having  a  conductivity  of  ",kp*241,"  Btu/hr-ft-F") ;} 
return (0); 

}  /% 

print() 

stream  *  fopen  ("LPT1",V) 
header (); 

fprintf  (stream, "\n%s" ,  "t,hr  Distance, in. ") ; 

fprintf  (stream,  "\n\r\017%s'\"  "}; 

for  (i=Q;i<*ST;i++) {fprintf  (stream,  "%2.1f\t",l*i*delx/2.54);} 
fprintf  (stream, "\n\r") ; 
for  (i*0;i<=8*(ST+2);i++) {fprintf  (stream, 
fprintf  (stream, "\n\r");  • 

switch (cres){ 
case  'a': 

for  (j=s0;j<s,nz;j++)  {fprintf 
(stream, "\n\r%2.2f  %s\t" ,  j*delt, " | ") ; 
for  ( i  =*0;  i<*ST ;  i  ++)  {fpri  ntf 
(stream, "%4. lf\t" , 1 .8* (tg*th[i] [j] -273)+32) ; }} 
break; 
case  *b': 

for  (j=0;j<=nz;j++) {fprintf 
(stream, "\n\r%2.2f  %s\t" , j*delt, " | ") ; 
for  (i=0;i<=$T;i++) {fprintf (stream, "%4.1f\t",cs[i] [j]);}} 
break; 
case  ' c 1 : 

for  (j*0,’j<nnz;j++){furintf 

(stream, "\n\r%2 . 2f  %s\t " ,  j*del t , " | " ) ; 

for  ( i =0 ; i <*ST ; i ++) { f p r i n t f ( s t ream , " %4 . 1 f \ t " , eg [ i ] [ j ] ) ; } } 

break; 

fprintf  (stream, "\r"); for  (i=l;i<=8*(ST+2);i++){ 

fprintf  (stream,  " 

fprintf  (stream, "\(522\014") ; 

fflush  (stream); 

return  )) ; 

^ard() 

I 

int  jk»0; 

while  (jk^lOP  0) {j k++; } 
return  (0) ; 


ovr() 

{ 
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} 


settextcolor(lO) ;  settextposition  (2,16); 

r  f  /rvr>i  i  4  \  r  _  . _ i/lin-xL  •  J _ 1 >  J.  - 


Tf  (BFL[ll=-ll){_outtext 
if  (BFL[2j==ll){j3uttext 
if  Bfl[3]==ll)-_outtext 
if  (BFL[4]==il) {outtext 
settextpositionl3,16) ; 

Tf  (BFL [8’’  . 

if  (BFL[9 
if  (BFLLOj 


"Both  Sides  Heated");' 
"One  Side  Air-Cooled"] 
"One  Side  Moist  Soil"! 
"One  Side  Insulated");}* 


=11) {_outtext(" Agent  =  HD");  • 

==11)  {  outtext ("Agent  =  GB"); '■ 

..  , _ .,.,==11) {“outtext ( " Agent  =  VX");! 

sprintf(bufl,"%3.1f  %s",coav,"  mg/g“); 
_settextposition(4,16);_outtext("Ave.  Init.  Conc.= 
return (0) ; 


11 )  ;_outtext(bufl) ; 


fnhp() 

switch (cres){ 
case  1  a' : 

hy= (1520+4309  'th[i][j]-255.3/tg)/(l-255.3/tg)); 
break; 
case  ' b ' : 

hy  =  1520+4309*(cs[i][j]/(4*coav)); 
break; 
case  ' c 1 : 

hy  =  1520+4309*cg[i][j]/(4*coav); 
break; 
case  1 d 1 : 
hy=0; 

for(i=Q;i<=ST-l; i++){ 

hy=hy+.5*(cs[i][j]+cs[i+l][j]+cg[ij[j]+cg[i+l]  [j])/ST;} 

hy=1520+4309*(hy/coav) ; 

break; 


prnscrnQ 


return (hy) ; 


stream  =  fopen("LPTl" ,"w") ;fprintf (stream, "\v\v\033\101\007") ; 
if  (BFL[5]==3){fprintf (stream, "\033\062") ;} 
for(j=0;j<=(int) (.914*ypix);j+=7){for(i=0;i<=639;i++){tm[il=0; 
for ( k=0 ; k<=7 ; k++) {m [k] =  getp i xel ( i , k+j ) ; i f (m[k]>l ) {m [k] =1 ; } 
tm[i]=tmri]+m[k3*pow(2,l7-k));}j 

fprintf  (stream, "\v\r\t\033\114\200\002"); for (i=0; i<=639; i++){ 
fputc(toascii (tm[i]) .stream) ;} 

fprintf (stream, "\v\v\r"); 

fprintf (stream, "\033\101\011") ; fprintf (stream, "\033\062") ; 
headerO ;fprintf (stream, "\014") ; 
fflush (stream) ; 
fcloseall; return (0) ; 


} 


